Vortex oscillations induced by a spin-polarized current in a magnetic nanopillar: 
Evidence for a failure of the Thiele approach 
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We investigate the vortex excitations induced by a spin-polarized current in a magnetic nanopillar 
by means of micromagnetic simulations and analytical calculations. Damped motion, stationary 
vortex rotation and the switching of the vortex core are successively observed for increasing values 
of the current. We demonstrate that even for small amplitude of the vortex motion, the analytical 
description based the classical Thiele approach can yield quantitatively and qualitatively unsound 
results. We suggest and validate a new analytical technique based on the calculation of the energy 
dissipation. 
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A magnetic vortex is a curling magnetization distribu- 
tion, with the magnetization pointing perpendicular to 
the plane within the nanometer size vortex core. This 
unique magnetic object has attracted much attention re- 
cently because of the fundamental interest to specific 
properties of such a nanoscale spin structure. Gyrotropic 
modes of vortices in magnetic nanocylinders have been 
intensively studied theoretically [l| and experimentally 
0. Apart from their fundamental relevance, the unique 
properties of the vortices are of considerable practical 
interest for new applications to memory and microwave 
technologies. In this view, the reversal of the magnetiza- 
tion within the vortex core by magnetic field and spin- 
polarized current has been thoroughly studied d, 0) H S] ■ 
More recently, sub- GHz dynamics of magnetic vortices 
induced by the spin transfer effect observed in nanopil- 
lars and nanocontacts 0, S, have raised a strong inter- 
est. Indeed, the associated microwave emissions in such 
vortex-based Spin- Transfer Nano-Oscillators (STNOs) 
occur at low current densities, without external magnetic 
field, together with high powers and narrow linewidths 
(< 1 MHz) comparatively to single-domain STOs. 

Traditionally, the analytical description of the vortex 
gyrotropic motion is based on the general approach for 
a translational motion of a magnetic soliton in an infi- 
nite media developed by A. Thiele [lo|. This calculation 
consists in a convolution of the Landau-Lifshitz-Gilbcrt 
(LLG) equation with the magnetization distribution un- 
der a specific condition of a translational motion of the 
magnetization pattern. Eventually a single equation (of- 
ten referred to as the Thiele equation) for the vortex core 
position X can be derived. The approach developed by 
Thiele to build his equation has been used for a long time 
to derive equations of vortex motion in many magnetic 
systems. In particular, it is often used to describe ana- 
lytically the vortex oscillations induced by spin current 
[llLll2Lll3l.ll4| in magnetic nanodiscs, in the ^ current per- 



micron discs can not be considered as translational due 
to a strong deformation of the vortex structure by the 
edges Guslienko et al. demonstrated that this de- 

formation should be taken into account in the calculation 
of the system energy [l| . However we show here that the 
impact of the spacial confinement on the vortex dynam- 
ics is much deeper. Our results evidence that, even if 
a proper model magnetization distribution is used, the 
Thiele approach applied for CPP nanodiscs with a spin 
current can give rise to significant qualitative and quan- 
titative errors. We suggest a new analytical technique to 
estimate the spin current-induced effects in such systems. 

In our calculations we consider vortex motion in a 
vortex-STNO. The system under study is sketched in the 
inset of Fig. [TJ The nanopillar spin valve has a circular 
cross-section. The reference layer is a fixed perpendicu- 
lar polarizer, the polarization vector p is perpendicular 
to the plane, and, to be clear in our interpretations, we 
disregard the stray magnetic field emitted by it. The 
initial magnetization distribution in the free layer is a 
vortex; the magnetization within the vortex core is par- 
allel to p. The current flow is assumed to be uniform 
in the pillar, with a radial symmetry of the current lines 
in the contact pads. The spin transfer term [16j in the 
calculations is given by (crJ/Afs)M x (M x p), J is the 
current density, M is the magnetization vector, is 
the magnetization of saturation and a represents the ef- 
ficiency of the spin-transfer torque: a = HP/ {2 |e| LMg), 
P is the spin polarization of the current, e is the charge 
of the electron, L is the sample thickness. 

As a starting point for our analytical calculations, we 
use a shortened form of the Thiele equation, which was 
used in Ref. 1] to account for the frequency of low- 
amplitude vortex oscillations in magnetic nanodiscs: 
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Here the gyrovector G is given by: 
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7 is the gyromagnetic ratio, 9, tp are the magnetization 
angles, and the integration in Eq. Q is over the magnetic 
disc. VK(X) is the potential energy of the shifted vortex. 
Guslienko et al. [l| showed that an appropriate model 
magnetization distribution for a moving vortex is given 
by the two-vortices ansatz (TVA): 
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Here p,x are polar coordinates in the disc plane, a = |X| 
is the vortex core displacement, (a, Xv) is the position 
of the vortex core center, R is the radius of the dot and 
C = 7r/2 or C = — 7r/2 for different regions of the dot. 
The TVA Eq. ([3]) defines a spin structure that satis- 
fies the magnetostatic boundary conditions, i.e. assumes 
zero magnetic charges on the side borders of the disc. 
The out-of-plane magnetization component Mg cos 9 can 
be described by a bell-shaped function, that is non-zero 
in the core region few nanometers in diameter [l^ . Using 
Eq. ([1]), hereafter wc address analytically the vortex gy- 
rotropic motion with a frequency f — Lo/in and a small 
amplitude a: 



X = 27r/e^ X X, a« R. 
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It follows from Eq. ([2]) that the gyrovector is given 
by G = — Gbz, where G = 2'kMsL/^ is the gyrocon- 
stant At small current densities, the major con- 

tribution to the vortex energy W(X) is the magneto- 
static energy Wm, arising from the appearance of vol- 
ume magnetic charges for a shifted vortex. It is given 

Recent simulations have 
shown that the contribution of the Oersted magnetic 
field generated by the current can be very important 
[3] ■ Therefore we also calculate the energy contribution 
Woe due to the Oersted field. Woe is given by integra- 
tion of the energy density — Hoe(r) • M(r,X) over the 
volume, where Hoe(r) is the Oersted field distribution 
at a given point r. The integration yields Woe (a) = 
l.TGnLRMsJa? / c [IS]. Summing up these two contribu- 
tions and using Eq. ([!]) and (|4]), one gets the analytical 
prediction for the vortex frequency: 



/ — /o" + foeJ, 



(5) 



where /q™ = §jM,L/R and foe = '-^jR/c. 

The critical current to launch the vortex oscilla- 
tions can be found considering the energy balance in 



the system. The energy dissipation is given hy W = 
/ (H ^ + If ' where || and || can be substituted by 
expressions found from the LLG equation as done in Ref. 

Assuming that outside of the core region 9 ^ n/2 
and neglecting the contribution from the core, one finds: 



W 



-a ip aV 

1 



M, 



iaJ)ipdV, (6) 



a is the Gilbert damping. The first term in Eq. ([6|) 
corresponds to the energy loss due to the damping, and 
the second term to the dissipation arising from the spin 
transfer torque, which can be positive or negative accord- 
ing to the current sign. For a steady motion the energy 
is conserved {W = 0), thus after some algebra, one can 
find from Eq. (g]): 
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where ?/ = ^ln{R/2ce) — |, Ce is an integration limit, that 
is of the order of exchange length 1^ = ^/ A/At:M^ [l^] , 
A is the exchange constant. From Eq. ([5]) and ^ we 
get an expression for the critical value Jci of the current 
density to excite the vortex oscillations: 
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A different prediction follows from the Thiele ap- 
proach. The Thiele equation that takes into account the 
damping and the spin transfer effect in the CPP config- 
uration is given by 111, 1^ : 
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where the spin transfer force Fst is: 
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(a different expression for Fst has been derived for GIF 
systems [l3|)- D is the damping tensor: 
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For circular dots, D = DE, where £^ is a unit tensor and 
the damping constant is D = arj'G [13]. The factor rj' 
and the previously introduced rj define the same quantity 
even if they are given by different expressions as we dis- 
cuss below. Calculation of the spin transfer force for the 
TVA yields Fst = 2nMsL(TJae^ [ill, El- For a steady 
gyrotropic motion, the third and the last terms of Eq. 
([9|) are perpendicular to the first and the second terms. 
Therefore the frequency of the vortex motion, given by 
Eq. ([T]), is not affected by the new terms of Eq. ([5]); 
instead, they define the amplitude of the vortex motion 
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llj. For the steady motion, the damping term is bal- 
anced by FsT, thus one gets: 
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Comparing this to Eq. ([7]), we see that Eq. (|12p yields 
about a twice smaller value of the first critical current 
Jci- The reason of this difference is related to the break- 
down of the assumption of a translational motion for the 
vortex, thus of the basic underlying assumption for the 
Thiele approach. 

Indeed, the derivation of Eq. 0, (fTT|) and PH)) essen- 
tially uses the following feature of the translational mo- 
tion of a magnetic soliton: M = — ^X, M 11 1. For a 

vortex moving in a nanodisc, however, the left-hand side 
of this expression vanishes at the disc side border, but 
the right-hand side has a finite value. It can be shown 
that for a << R, most contribution to the magnitude 
of the gyroforce G and the damping D comes from the 
vicinity of the disc center (p < R/2); the magnetization 
motion within his region can be considered translational, 
to some approximation. In contrast to it, even for a very 
small a, about a half of the magnitude of Fst given by 
Eq. (fTO|) originates from the region of the disc close to 
its boundary. 

Our analytical results can be compared to numerical 
micromagnetic simulations. In the simulations a nanopil- 
lar 300 nm in diameter is considered. The free layer is 10 
nm thick and has the following magnetic parameters: Ms 
= 800 emu/cm^, A = 1.3 x lO"*" erg/cm and a = 0.01 
(values for permalloy). We use a two-dimensional mesh 
with in-plane cell size 1.5 x 1.5 nm^. The polarization is 
taken to be P = 0.2. The micromagnetic simulations are 
performed by numerical integration of the LLG equation 
using our micromagnetic code based on the forth order 
Runge-Kutta method with an adaptive time-step control 
for the time integration. 

We observe vortex excitations only for positive current, 
which is defined as a flow of electrons from the free layer 
to the polarizer. The vortex motion is damped for small 
current densities J < Jci, where the first critical current 
density Jci = 4.9 x 10^ A/cm^. For larger currents, 
after some transitional period, the vortex is gyrating on 
a steady circular orbit. Interestingly, Jci is about one 
order of magnitude less than the critical current density 
for excitation of magnetization oscillations in nanopil- 
lar STNOs with nominally uniform free magnetic layer 



211 . 122] . The values of J, for which the steady vortex os- 
cillations are observed, are limited by the second critical 
current value Jc2 =9.0 x 10^ A/cm^. For J > Jc2, on 
reaching a critical orbit, the core of the vortex is reversed. 
The details of this process: appearance of a vortex with 
opposite polarity and an antivortex, annihilation of the 
latter with the original vortex, essentially reproduce the 
previous finding for the vortex core switching by the 
field or current [3, S, [i2| ■ After the reversal of the core. 



the direction of the vortex gyration is changed and the 
vortex oscillations are damped. 

For each point within Jci < J < Jc2, the vortex mo- 
tion is simulated for 100 ns after reaching a stationary 
orbit. The vortex frequencies extracted from these sim- 
ulations together with the radius of the oscillation orbit 
are presented in Fig[TJ On increasing J, the oscillation 
frequency increases ranging from 0.34 GHz to 0.41 GHz. 
The radius of the orbit increases with the current as well, 
reaching 125 nm at J = 8.5 x 10^ A/cm^. 
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FIG. 1: (color online) Steady vortex gyration induced by the 
spin polarized current: frequency / = ui/2tt (squares) and 
radius of the vortex core orbit a (triangles) as a function of 
the current density J (numerical simulations). Top insets: 
illustrations for the time evolution of the averaged projection 
of the free layer magnetization on the polar axis, for J < 
Jci (left), Jci < J < Jc2 (center) and J > Jc2 (right), in 
arbitrary units. Bottom inset: sketch of the device geometry. 

Analytical prediction for the vortex frequency at J = 
Jci, given by Eq. ([5|), is f = 0.36 GHz that is in a good 
correspondance to the simulation results / = 0.34 GHz. 
This is in agreement with the previous results of Ref. 
[l|, which showed that Eq. ([T]) can properly describe 
the frequency of small-amplitude vortex oscillations in 
magnetic nanodots. The factor rj for our system can be 
extracted using additional simulations. We find that if 
the current is switched off, the vortex motion is a gyration 
with the orbit damped in time as a oc exp(— t/r), r is 
a time constant. On the other hand, it follows from Eg . 
([6]) and (jl]) that at zero current a cx exp {—2aTTrift) [23|. 
Comparing this to the numerical results, we find rj = rj' 
= 1.65. This value is in reasonable agreement with the 
prediction of the current work (ry = 1.34 taking Cf, = le) 
and previously made analytical estimations (e.g. ry' = 2.3 
from Ref. fl). 

Taking 77 and / from the simulations, we find from 
Eq. jl]) that the analytical prediction Jci = 4.5 x 10^ 
A/cm^ is in good agreement with the numerical result. 
The resuh of the Thiele approach (Eq. ^) is Jci = 2.3 
X 10^ A/cni^, that illustrates our statement of its failure 
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to account for the vortex motion in CPP nanopillars. 
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FIG. 2: (color online) Numerical result for the constrained 
polarizer, in the notations of Fig. [1] Inset: sketch of the 
device geometry. 

We demonstrate the breakdown of the Thiele approach 
in a different way. We perform simulations for a config- 
uration that we call constrained polarizer, for which the 
current polarization P equals to 0.2 for p < 50 nm and 
P = 0, hence cr = 0, for larger p. Thus in contrast to 
the case of the uniform polarizer for this new structure, 
the current does not excite the regions close to the disc 
boundary. The resulted dependences /(J) and a(J) are 
presented in Fig. [21 We find for this configuration that 
J(ji = (4.9 ± 0.1) xlO^ A/cm^ equals to the critical cur- 
rent for a uniform polarizer. This is in perfect agreement 
with the analytical result from the energy consideration. 
Indeed, Eq. ^ gives equal results, hence equal values of 
Jci, for both configurations aX a « R. 

The prediction of the Thiele approach for a constrained 
polarizer is that the magnitude of the spin-transfer force 
FsT (Eq- pO)) ) is a factor 1.8 less than that for a uniform 
polarizer. Accordingly, the value of Jci in the two config- 
urations differs at about the same factor, in contradiction 
to the numerical results. Thus we demonstrate that the 
Thiele approach can give rise to not only quantitative, 
but also qualitative disagreements. 

Simulation results for the constrained polarizer contain 
other remarkable facts. We see that the second critical 
current Jc2 = 6.5 xlO'^ A/cw? is by a factor of about 7 
larger that Jc2 for a uniform polarizer. The oscillation 
orbit gradually increases with the current reaching a = 
108 nm at J = Jc2- The oscillation frequency starts 
at f = 0.34 GHz at Jci like for a uniform polarizer but 
reaches a larger frequency f = 0.52 GHz at J = Jc2- 

It has been predicted that the vortex core is reversed if 
its velocity reaches a critical value Vcrit that is about 340 
m/s for permalloy independently of the device design 
Our simulation results are 320 m/s for the uniform po- 
larizer and 350 m/s for the constrained polarizer in nice 



agreement with this prediction. A small difference be- 
tween the values is presumably due to a different extent 
of the vortex deformation at critical orbits. The value 
of Vcrit together with the dependencies of /(J) and a( J) 
define the value of the second critical current Jc2- As for 
the constrained polarizer /(J) and a(J) are much less 
steep functions than those for the uniform polarizer, Jc2 
in the former is substantially larger than in the latter. 
Larger frequencies at a given orbit, found for the con- 
strained polarizer, are related to the strong influence of 
the Oersted field. These facts make this configuration 
promising for applications. It can be implemented by 
reducing the polarizer dimensions or simply by using a 
point contact technique. 

In conclusion, we clearly show that the Thiele formal- 
ism fails to give a sound analytical description for the 
vortex motion in the CPP magnetic nanopillars due to 
the breakdown of the underlying assumption of the trans- 
lational motion of the vortex. We suggest a new original 
analytical approach based on calculation of the energy 
dissipation, which is shown to be in good agreement with 
the numerical results. Our calculations demonstrate that 
vortex-based STNOs can potentially have very low values 
of J(7i, large values of Jc2 and operate at zero external 
magnetic field, which make them promising candidates 
for future microwave technology applications. 
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